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ABSTRACT 

The eclipsing novalike cataclysmic variable star V348 Pup exhibits a persistent lumi- 
nosity modulation with a period 6 per cent longer than its 2.44 hour orbital-period 
(-Porb)- This has been interpreted as a "positive super hump" resulting from a slowly 
precessing non-axisymmetric accretion disc gravitationally interacting with the sec- 
ondary. We find a clear modulation of mid-eclipse times on the superhump period, 
which agrees well with the predictions of a simple precessing eccentric disc model. 
Our modelling shows that the disc light centre is on the far side of the disc from the 
donor star when the superhump reaches maximum light. This phasing suggests a link 
between superhumps in V348 Pup and late superhumps in SU UMa systems. Mod- 
elling of the full lightcurve and maximum entropy eclipse mapping both show that 
the disc emission is concentrated closer to the white dwarf at superhump maximum 
than at superhump minimum. We detect additional signals consistent with the beat 
periods between the implied disc precession period and both ^Porb and |P rb- 

Key words: accretion discs - binaries: close - binaries: eclipsing - stars: individual: 
V348 Pup - stars: cataclysmic variables 



1 INTRODUCTION 

V348 Pup (1H 0709-360, Pup 1) is a novalike cataclysmic 
variable (CV): a system with a high mass transfer rate which 
maintains its accretion disc in the hot, ionized, high viscosity 
state reached by dwarf novae in outburst. It exhibits deep 
eclipses in its optical and infrared lightcurves (Tuohy et al. 
1990): it is a high inclination system with orbital period 
Porb = 2.44 hours (Baptista et al. 1996). 



1.1 Superhumps 

Modulations in luminosity with a period a few per cent 
longer than the orbital period have been observed in many 
short period CVs (see reviews in Molnar & Kobulnicky 1992, 
Warner 1995, Patterson 1998a). These modulations typically 
take the form of a distinct increase in luminosity, or super- 
hump. The standard explanation of this phenonemon is that 
the system contains an eccentric precessing accretion disc. 

If the accretion disc extends out far enough, the out- 
ermost orbits of disc matter can resonate with the tidal 
influence of the secondary star as it orbits the system. A 
3:1 resonance can occur which results in the disc becom- 
ing distorted to form an eccentric non-axisymmetric shape. 
The tidal forces acting on this eccentric disc will cause it to 
precess slowly in a prograde direction. 

The superhump period, P s h, is then the beat period 



period, P or b (Osaki 1996): 



1 
PLh" 



P a h is the period on which the relative orientation of the line 
of centres of the two stars and the eccentric disc repeats. Pos- 
sible models for the light modulation on P s h are described 
below. This paper considers these models in relation to our 
observations. 

In the tidal model the superhump is a result of tidal 
stresses acting on the precessing eccentric disc (Whitehurst 
1998b). The light may be due to a perturbation of the veloc- 
ity field in the outer disc, leading to azimuthal velocity gra- 
dients and crossing or converging particle trajectories. Thus 
extra dissipation modulated on the superhump period arises 
when the secondary sweeps past the eccentric disc. In ad- 
dition, the superhump-modulated tidal stress would lead to 
a superhump-modulated angular momentum loss from the 
disc which would facilitate a variation in the mass trans- 
fer rate through the disc, and hence a modulation in disc 
luminosity. 

The bright spot model arises from noting that the en- 
ergy gained by material in the accretion stream will depend 
on how far it falls before impacting on the disc (Vogt 1981). 
The energy dissipated at impact will be modulated on the 
superhump period since the non-axisymmetric disc radius 
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causes a stream-disc impact region at varying depths in the 
white dwarf potential well. 

Recent SPH simulations of accretions discs in AM CVn 
stars lead to a third, more realistic, model in which the disc 
shape changes from nearly circular to highly eccentric over 
the course of a superhump period (Simpson & Wood 1998). 
Superhumps arise from viscous energy production as the dis- 
torting disc is tidally stressed. Other SPH simulations (e.g. 
Murray 1996, 1998) also reveal a disc whose shape changes, 
with Murray (1996) predicting superhump modulations from 
both the periodic compression of the eccentric disc and the 
varying depth in the primary Roche potential at which the 
stream impacts the disc. 

Dwarf novae in super-outburst exhibit two distinct pos- 
itive superhump phenomena (Vogt 1983, Schoembs 1986). 
Normal superhumps appear early in the super-outburst and 
fade away towards the end of the outburst plateau to be 
replaced with 'late' superhumps which persist into quies- 
cence. These late superhumps are roughly anti-phased with 
the normal superhumps, and are more likely to be analogous 
to the persistent superhumps seen in novalikes (Patterson 
1998b), where the system has had sufficient time to settle 
into a steady state. Our extensive photometry (Section ^ 
reveals similarities between superhumps in V348 Pup and 
late superhump s in dwarf novae. 
In Section 
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we present power spectra revealing the 
superhump period and addi tion al signals close to orbital pe- 
riod harmonics. In Section 3.2 we estimate the orbital pa- 
rameters, q and i, for V348 Pup using the average orbital 
lightcurve and the superhump period. The wav efor m of the 
sup erhump modulation is discussed in Section 3.3. Section 
3.4 considers average orbital lightcurves grouped according 
to superhump phase. In Section we fit our lightcurves 
with a precessing eccentric disc model, hence deducing the 
location of light centre of the disc. We consider the results 
of maximum entropy eclipse mapping in section 3.6. Our 



results and their implications are discussed in Section 



2 THE OBSERVATIONS 

The observing campaign comprises 24 nights of rapid pho- 
tometry from December 1991, February 1993 and January 
1995 (see Table 0). The 1991 and 1993 observations (12 and 
8 nights respectively) were taken using the the 40-inch tele- 
scope at CTIO with a blue copper sulphate filter. The Jan- 
uary 1995 run consists of 4 nights of R band data. All the 
data have been corrected for atmospheric extinction and the 
1995 data have also been calibrated to give an absolute flux. 
In 1995 the average out of eclipse R magnitude is 15.5 mag; 
at mid-eclipse R = 16.8. Examples of typical data are plot- 
ted in Figure hi. 



3 ANALYSIS 

We determined mid-eclipse timings as described in Section 
3.5 from which we determined an orbital ephemeris for this 
dataset 

T mid = HJD 2448591.667969(85) 



+ 0.101838931(14)E. 
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Figure 1. 4 cycles of typical data from 1991 and 1993. 



Table 1. Log of observations 



Date 


HJD start 


Duration 


Integration 


Filter 




-2440000 


(hours) 


time (sec) 




1991 Dec 1 


8591.6779 


3.96 


42 


CuS0 4 


1991 Dec 2 


8592.6138 


3.55 


42 


CuS0 4 


1991 Dec 3 


8593.6092 


5.68 


42 


CuS0 4 


1991 Dec 4 


8594.6249 


5.14 


42 


CuS0 4 


1991 Dec 5 


8595.6203 


5.37 


42 


CuS0 4 


1991 Dec 6 


8596.6328 


5.18 


42 


CuS0 4 


1991 Dec 7 


8597.6293 


5.08 


42 


CuS0 4 


1991 Dec 8 


8598.6107 


5.59 


42 


CuS0 4 


1993 Feb 2 


9020.5672 


5.56 


10 


CuS0 4 


1993 Feb 3 


9021.5549 


6.24 


10 


CuS0 4 


1993 Feb 4 


9022.5497 


6.17 


10 


CuS0 4 


1993 Feb 9 


9027.5565 


3.71 


10 


CuS0 4 


1993 Feb 11 


9029.7051 


2.36 


10 


CuS0 4 


1993 Feb 12 


9030.6851 


2.84 


10 


CuS0 4 


1993 Feb 13 


9031.5446 


6.02 


10 


CuS0 4 


1993 Feb 14 


9032.5427 


6.23 


10 


CuS0 4 


1993 Feb 21 


9039.6498 


3.20 


10 


CuS0 4 


1993 Feb 22 


9040.5572 


3.13 


10 


CuS0 4 


1993 Feb 23 


9041.5479 


2.52 


10 


CuS0 4 


1995 Jan 3 


9720.5689 


4.83 


10 


R band 


1995 Jan 4 


9721.5521 


3.32 


10 


R band 


1995 Jan 5 


9722.5509 


7.15 


10 


R band 


1995 Jan 6 


9723.5472 


7.24 


10 


R band 



This is consistent with the Baptista et al. (1996) ephemeris 
within the quoted error limits. We adopt our ephemeris for 
this analysis. The eclipse timings are given in Table 0. 



3.1 The superhump period 

Before performing any analysis, we normalized each night of 
data by dividing by the average out of eclipse value. 

To make detection of non-orbital modulations in the 
data easier, the average orbital lightcurves^] from each years' 
observations were calculated and subtracted from the corre- 
sponding data. The resulting lightcurves contain no orbital 

* Superhump phase grouped average lightcurves are shown in 
Figure Pj. 
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Table 2. Eclipse timings 

HJD mid-eclipse Error 
-2440000 estimate 



HJD mid-eclipse 
-2440000 



Error 
estimate 



December 1991 




February 1993 


8591 


769462 





000121 


9020 


613001 





000036 


8592 


686379 





000096 


9020 


715111 





000043 


8593 


705382 





000120 


9021 


632451 





000090 


8593 


806642 





000017 


9021 


734218 





000082 


8594 


722715 





000090 


9022 


650034 





000232 


8594 


824510 





000162 


9022 


752227 





000175 


8595 


640158 





000104 


9027 


640565 





000417 


8595 


742287 





000201 


9029 


779356 





000091 


8596 


657488 





000082 


9030 


695522 





000026 


8596 


759320 





000100 


9030 


797304 





000213 


8597 


676629 





000068 


9031 


612099 





000105 


8597 


778677 





000083 


9031 


714122 





000104 


8598 


694383 





000058 


9032 


630438 





000045 


8598 


796652 





000090 


9032 


732446 





000022 




January 


1995 




9039 


657734 





000007 


9720 


653953 





000003 


9039 


759941 





000159 


9720 


756257 





000300 


9040 


573722 





000065 


9721 


571129 





000069 


9040 


675596 





000025 


9721 


672783 





000132 


9041 


592531 





000082 


9722 


589171 





000043 










9722 


690736 





000053 










9723 


607919 





000087 










9723 


709706 





000135 










9723 


811139 





000068 











variations. Lomb-Scargle periodograms were calculated for 
each year's data, and are shown in Figures |and|[ The 1991 
periodogram reveals a periodicity with period 0.10763 days 
and simple aliasing structure. The 1993 periodogram has 
higher resolution, more complicated alias structure, and the 
strongest periodicity at 0.10857 days. The 1995 power spec- 
trum (having the least time coverage) is lower resolution, 
however a clear signal at 0.10760 days and its aliases is 
present. By comparison with the clearer 1991 and 1993 spec- 
tra we surmise the 0.10760-day peak is the true signal. These 
periods are all close to 6 per cent greater than P or b and cor- 
respond to the period excesses, e, shown in Table |j| defined 
using superhump period, P a h = (l + e)Port>- The inferred disc 
precession periods are also shown in Table |§| 

It is notoriously difficult to determine errors in peri- 
ods measured from periodograms. To estimate the errors in 
the periods detected here, various fake datasets were gener- 
ated. The lightcurves were smoothed and the residuals used 
to characterize the variance of the random noise. Various 
types of noise with the same variance were added to both 
the smoothed lightcurves and smoothed superhump modu- 
lations. The variance in the period value determined from 
these different datasets provides a measure of the error in 
the period measurement. Since the smoothed data will still 
contain noise artefacts these are probably underestimates. 
The 1993 dataset with its 20-day time base should provide 
the most precise period. Therefore, the possibility that the 
real 1991 and 1995 superhump periods are in fact closer to 
the 1993 value than our measurements suggest should not be 
ruled out, as Figure [] shows. However, variability in the de- 
tected periodicity has also been seen before for superhumps 
(Patterson 1998b). 

Superhumping systems show a clear correlation between 




Figure 2. Power spectra. Arrows from left to right indicate the 
superhump period and the sidebands of the first and second har- 
monics of the orbital period (see Table U). 



e and P OT b, with e increasing with P or t>; the superhump pe- 
riods we detect in V348 Pup are consistent with this trend 
(Figure 3, Patterson 1998b). 

The power spectra also reveal signals at frequencies cor- 
responding to sidebands of harmonics of the orbital period. 
The strongest such detections are at periods corresponding 
to 2f2 or ;, — Qprec and 3f2 or 5 — f2 pre c (marked with arrows in 
Figure ^|): the predicted and the directly measured values 
of these sidebands are shown in Table ^, and graphically in 
Figure ^. The errors in these periods were estimated using 
the method described earlier. The detected periods do not 
all agree to within the estimated errors, suggesting that the 
error estimates may be a little too low, as expected; the 
highest quality 1993 data agrees best. 

The simplest way to produce these sidebands is by mod- 
ulating the brightness or visibility of the superhump with 
orbital phase. If we consider the orbital lightcurve as a sum 
of Fourier components with frequencies n£l or b, then follow- 
ing the approach of Warner (1986) and Norton, Beardmore 
& Taylor (1996), the eclipse of the superhump light source 
will produce signals at frequencies (n + l)Q, or b — il pr ec- Sig- 
nals at frequencies (n — l)fi or t, + Q pre c are also predicted 
but we find no evidence of these; perhaps they are nulli- 
fied by other signals of the same frequency in antiphase. In 
Section 3.4 we present evidence for a correlation between 



superhump amplitude and orbital phase; this modulation of 
superhump amplitude with orbital phase could also lead to 
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December 1991 



February 1993 




Table 4. Sidebands of the orbital period harmonics 



Frequency (cycles/day) 



Figure 3. Detail of the power spectra shown in Figure ^| Arrows 
indicate the superhump period. The orbital period P or b is also 
marked, from which it is clear that the orbital variation has been 
successfully removed. 



Table 3. Measured superhump characteristics, period P s /,, the 
period excesses e, the inferred disc precession periods P prec and 
the consequent estimates of the mass ratio q. 



Psh (d) 



(d) 



Fractional 
amplitude 



1991 0.107628(8) 0.0568 1.89 0.20 0.26 

1993 0.108567(2) 0.0661 1.64 0.15 0.31 

1995 0.10760(7) 0.0566 1.90 0.13 0.26 



the observed sideband signals. The SPH models of Simp- 
son & Wood (1998) predict the formation of double armed 
spiral density waves in the disc whose rotation rate, they 
suggest, might lead to observed signals at about three times 
the superhump frequency. They also suggest that viewing 
these structures from non-zero inclination could lead to the 
detection of further frequencies, although they do not make 
precise predictions. Observations of the dwarf nova IP Peg 
in outburst have revealed evidence of such spiral structure 
(Steeghs, Harlaftis & Home 1997). 

There is no significant signal around period P prec =1.6 
- 1.9 days in either the normalized or un-normalized 
lightcurves. The disc precession period is similarly absent 
in other persistently superhumping systems, notably AM 
CVn where P pre c is clearly revealed by absorption line spec- 
troscopy (Patterson, Halpern and Shambrook 1993). 



Pred. P (d) 



Qprec 

Meas. P 



(d) Pred. P (d) Meas. P (d) 



91 0.052327(1) 
93 0.0525477(4) 
95 0.05232(1) 



0.05241(1) 
0.052540(1) 
0.052493(8) 



0.0345661(6) 
0.0346623(1) 
0.034563(6) 



0.034617(1) 
0.03466(1) 
0.034588(2) 



19.03 19.06 19.09 19.12 

Frequency (cycles /day) 



28.83 28.88 28.92 28.5 
Frequency (cycles/day) 



Figure 4. Sidebands of the orbital period harmonics. Dotted lines 
mark predicted frequency range for each signal, solid lines show 
the measured range. 



3.2 The orbital parameters 

The mass ratio, q, of a superhumper can be estimated, given 
its period excess, e (Patterson 1998b). 

_ 0.23g 
£ ~ l+0.27g' 

This leads to the estimates of q for V348 Pup shown in Table 
^[ We favour the mass ratio, q = 0.31, estimated from the 
most accurate 1993 superhump period. We note that SPH 
simulations of eccentric discs (Murray 1998, 1999) suggest a 
more complicated relationship between e and q with the disc 
precession rate depending on the gas pressure and viscosity 
of the disc in addition to the mass ratio. This does not affect 
our substantive results. 

If we assume that the secondary star is Roche lobe fill- 
ing, the width, w, of eclipse of a point source at the centre 
of the compact object uniquely defines orbital inclination, i, 
as a function of q. We can thus compute i as a function of q 
and w. 

When the centre of the compact object (point P) is 
first eclipsed (orbital phase cii), about half of the disc area 
will be eclipsed, and therefore for a disc whose intensity 
distribution is symmetric about the line of centres of the 
two stars, the fraction of disc flux eclipsed at this phase will 
be ~ 0.5 (Figure ^). Similarly at the end of the eclipse of P 
(orbital phase (j>2) the fraction of disc light visible is again 
~ 0.5. Further assuming that the lightcurve consists purely 
of emission from a disc in the orbital plane, the full width 
of eclipse at half intensity will be equal to w. 

Using the average eclipses from each of our datasets 
to give w, we obtained i. We checked the assumption that 
half of the disc area is eclipsed at <j>\ and 4>2- for these values 
obtained it is a good assumption. We therefore adopt orbital 
parameters g = 0.31, i = 81? 1 ± 1° The conclusions drawn 
later are identical to those obtained using q = 0.36 with 
corresponding i = 80?0 ± 1° which resulted from an earlier 
estimate of q. 
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Table 5. Groups 

December 1991 February 1993 January 1995 





Pprcc 


= 1.89d 


Pprcc 


= 1.64d 


-fprcc 


= 1.90<f 


Group 


4>sb 




</>sh 




0sh 




A 


0.24 


0.06 


0.02 


0.07 


0.03 


0.06 


B 


0.71 


0.07 


0.48 


0.11 


0.54 


0.04 




Figure 5. The eclipsed region of the orbital plane at phase <j>l 
showing how approximately half of the area of a 0.9Rl disc is 
eclipsed at this phase. 



3.3 The superhump modulation 

To define a zero point in superhump phase each set of data 
was folded onto its detected superhump period, binned into 
100 phase bins, and a sine- wave was fitted to each (see Figure 

We assessed the contribution of flickering to these 
curves by using various methods of binning the data. Fig- 
ure ^| shows the curves produced from simply averaging the 
points in each bin. Since flickering in the lightcurve consists 
of brief increases in luminosity, by giving more weight to 
the lower values in each bin or using only the lower points 
in a bin, the impact of flickering on these superhump phase 
binned lightcurves should be reduced. We therefore gener- 
ated curves by averaging only the lowest 25 per cent of fluxes 
in each bin, and by weighting the lower values more strongly 
than the higher values. The curves produced by these differ- 
ent methods are virtually identical, except for a flux offset 
between each. This suggests that flickering has little effect 
on the superhump curves. Since we have extensive datasets 
and the timescales of the flickering and the superhumps are 
very different, this is not unexpected. 

For a high inclination system a modulation on the su- 
perhump period will arise due to the eclipses of the process- 
ing accretion disc changing as the disc orientation changes. 
This effect will occur in addition to the intrinsic variations 
in luminosity which are observed in non-eclipsing systems. 
Superhump modulations calculated using only non-eclipsing 
phases are almost identical to those in Figure ^, which im- 
plies that in this system the form of the superhump light 
curve is not affected by the changing eclipse shape. 

The broad form of the modulation is consistent for all 
three sets of observations. The peak-to-peak fractional am- 
plitudes are shown in Table |§| and decline steadily from year 
to year. The more detailed structure, particularly the re- 
gion with lower flux between (j> s h ~ 0.8 and 1.0 in the 1991 
modulation, appears to be genuine; we find neither flickering 
nor eclipses have significant effect. Furthermore the abrupt 
changes do not correspond to changes in system brightness 
from one night to the next. 

Simpson & Wood (1998) calculate pseudo- lightcurves 
for superhumps. Assuming the light emitted from the disc 
is proportional to changes in the total internal energy of 
the gas in the disc, they present superhump shapes calcu- 
lated for mass ratios of 0.050, 0.075 and 0.100 (their Figure 
5). These curves have significant differences in morphology: 



the q — 0.050 curve has a sharp rise and slow decline, the 
q — 0.075 curve is reasonably symmetric, and the q — 0.100 
curve has a slow rise and steeper decline. The cleanest and 
most reliable of our superhump curves, that from 1993, also 
shows an asymmetric shape, with a slow rise and sharper 
decline, agreeing best with their highest value of q = 0.100. 
We expect q ~ 0.31 so this is encouraging. 



3.4 Average lightcurves 

Because of the sampling of our data (Table [j]), the observed 
eclipses for each year appear in two rough groups, with su- 
perhump phases separated by about 0.5. We compared the 
average orbital lightcurves corresponding to each group. The 
average mid-eclipse superhump phase, </> 8 h, of each group, 
and the range as indicated by the variance of </> a h for each 
group are shown in Table ^, the average orbital curves are 
shown in Figure ^| 

The group A curve for 1991 (Figure displays a clear 
hump at orbital phase around -0.2: the mid-eclipse super- 
hump phase of group A is 0.24, so eclipse should occur 0.25 
in orbital phase after a superhump maximum. Group B has 
mid-eclipse superhump phase 0.71, meaning that the eclipse 
should occur around 0.3 in orbital phase before the super- 
hump. The post-eclipse flux in group B is higher than in A, 
though the hump is not so sharp as in A. 

Group A in the 1993 data displays a hump peaking just 
before mid-eclipse, while the curve for group B is rather flat 
out of eclipse. The mid-eclipse superhump phase of group 
A is 0.02, while group B has superhump phase 0.48. This 
is again consistent with the position of the superhump. The 
difference curve (i.e. A-B) seems to show a broad super- 
hump which is partially eclipsed by the secondary. We ex- 
pect group B to display a hump around orbital phase 0.5. 
This is not obvious, meaning that the superhump is more 
prominent when is occurs at orbital phase (group A) than 
at phase 0.5. The superhump light is therefore not emitted 
isotropically£j Schoembs (1986) also noted a similar effect 
in OY Car. This will be considered further in Section 

The superhump phasing of the 1995 eclipses is almost 
the same as those observed in 1993, but the smaller extent of 
the 1995 data and the low fractional amplitude makes identi- 
fying the hump difficult without first subtracting the orbital 
light curve (compare Figs ^b and ^|c). However, the flux dur- 
ing eclipse for group A is higher than for group B, consistent 
with the phasing which suggests that a superhump should 
occur at mid-eclipse. 



t Eclipses of the superhump light will obviously be most impor- 
tant when the secondary is at inferior conjunction, i.e. when the 
superhump maximum is at orbital phase 
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December 1991 February 1993 January 1995 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Superhump phase, 3h Superhump phase, 3h Superhump phase, 3h 



Figure 6. The superhump modulations. For each dataset, the average orbital light curve was subtracted and the resulting data folded 
and binned on superhump phase. Dotted lines show the sine waves fitted to determine phasing 



3.5 Eclipse parameters 

The O-C mid-eclipse times and the eclipse widths are shown 
in Figure The mid-eclipse times were determined both by 
fitting a parabola to the deepest half of each eclipse and also 
by finding the centroid. The discrepancy between the two 
determinations provides an indication of the uncertainty. 

As the eccentric disc precesses slowly in our frame, we 
expect to see the eclipse width and midpoint phase mod- 
ulated on the apsidal precession period. These quantities 
will be similarly modulated in superhump phase, since the 
superhump phase and precession phase of an eclipse are 
both measures of the relative orientation of disc and sec- 
ondary star at mid-eclipse. Figure ^ shows that eclipse tim- 
ings for all years exhibit a precession period modulation. 
The widths also show evidence of a modulation. The limited 
superhump phase coverage means that conclusions cannot 
easily be drawn from inspecting the datapoints alone. Such 
variations in eclipse asymmetry have been observed in other 
superhumping systems e.g. OY Car (Schoembs 1986) and Z 
Cha (Warner and O'Donoghue 1988, Kuulkers et al. 1991). 
To further investigate the disc shape we produced a simple 
model which was then fitted to the observed lightcurves. 

Our simple eccentric disc prescription has a circular in- 
ner boundary with radius r m i„, centred on the white dwarf. 
The outer boundary is an ellipse of semi-major axis a m ax, ec- 
centricity e, with one focus also centred on the white dwarf. 
The disc brightness at distance r(a) from the white dwarf 
at an angle a to the semi-major axis is 



5(a) oc 



' ma x 

(a) - 

1 77 



where r maa; (a) is the distance from the white dwarf to 
the outer disc boundary at orientation a. Brightness con- 
tours are therefore circular at the inner boundary, smoothly 
changing to elliptical at the outer boundary. This form for 
S(a) reduces to 5 oc r~ n if the disc is circular. Our model is 
sensible for a tidally distorted disc, since the tidal influence 
of the secondary star is unimportant at the inner disc, so we 
expect a more or less circular inner disc. 

In an inertial frame, the disc slowly precesses progradely 
with period P pre c- With respect to the corotating frame of 
the system, this disc then rotates retrogradely with period 
P s h- Let the relative orientation of the line of apsides of the 
disc with the line of centres of the binary when superhump 
maximum occurs be (fidisc- The structure of the disc in our 



model is therefore described by five parameters : r m i n , a max , 
e, n and (j>disc- This model was chosen as the simplest way to 
model an eccentric precessing disc, with as few parameters 
as possible. A similar model was used by Patterson, Halpern 
& Shambrook (1993) to model the disc in AM CVn. 

We generated synthetic lightcurves for eclipses of our 



model disc using the orbital parameters from Section 3.2 



By varying the parameters of our model to minimize the 
reduced Xr of the fit, we obtained a best fit of our model to 
our lightcurves. The smoothed superhump was subtracted 
from the lightcurves before fitting the model in order to 
remove the intrinsic variation in the disc flux, enabling us 
to study the shape of the disc. In Section B.4 we noted that 
the superhump is more visible for cf> s h = at 4> OT b = 0, so 
ideally we should subtract a superhump modulation which 
takes account of this variation in superhump prominence, 
but insufficient sampling of the disc precession phase by our 
data prevents us from doing this. 

The downhill simplex method for minimizing multidi- 
mensional functions was used (the AMOEBA routine from 
Press et al., 1996). Each orbit (centred on an eclipse) was 
allowed to have a different total disc flux and a different 
uneclipsed flux. This prevents the variation of ~10 per cent 
in flux from one orbit to the next from interfering with the 
results, and allows for the possibility of a contribution to the 
lightcurve from regions never eclipsed by the secondary star. 
The errors in the fluxes were estimated as being cr\J flux 
where a is the square root of the variance of the flux between 
orbital phase 0.2 and 0.8. This estimate therefore includes 
the effect of flickering. 

We used two methods to assess the robustness and ac- 
curacy of these fits. Monte-Carlo methods estimated the size 
of the region in parameter space which has a ~75 per cent 
chance of containing the 'true' values of the parameters. To 
assess the uniqueness and robustness of each solution we 
carried out the fitting process 20 times for each model/data 
combination starting each fit with a different random initial 
simplex. The solution chosen is that with the lowest value of 
reduced x ■ Extreme outlier solutions are rejected and the 
variance in the parameters for the remaining solutions is a 
measure of the accuracy with which the AMOEBA routine 
converges to a unique solution. The errors quoted in Table [?] 
for each parameter are whichever is the greater of confidence 
region estimate or the variance in the parameter from the 
multiple fits. The parameters resulting from all the fits are 
shown in Table [t| Rl is the Eggleton radius of the primary 
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Figure 8. Eclipse parameters. The O-C mid-eclipse times and the full eclipse widths at half depth. The superhump phase, </> s h, is 
calculated using the superhump period, P s h determined directly from the corresponding dataset. The continuous, da shed, dot-dashed 
and dot-dot-dot-dashed curves correspond to predictions from the lightcurve fits /, e, a and b described in Section |j.5| . 



Roche lobe (Eggleton 1983). The \ 2 surface in parameter 
space is not perfectly smooth. There is a broad global mini- 
mum superimposed with smaller amplitude bumps. Close to 
the global minimum the gradient of \ 2 is l° w an d so small 
bumps can lead AMOEBA to settle into a local minimum near 
the real minimum. 

Fits using the model as described above will be referred 
to as fit /. The emission extends out to 80 - 90 per cent 
of Roche lobe size, while the largest radius of the disc is 
dmaxi^- + e) = 0.97-Rl for 1995 data. These results suggest 
that the disc does indeed extend out to the tidal cut-off 
radius, r t id e ; r U de ~ 0.9Rl (Paczynski 1977). 

The 1995 radius is 10 per cent larger than the 1991 
and 1993 radii. The 1991 and 1993 observations are in blue 
CuSCU filter light, which should come from hotter inner disc 
regions; the 1995 observations are in R which is expected to 
weight the outer disc emission more heavily, perhaps causing 
the inferred disc radius to be largest for the 1995 observa- 
tions. A simple model black-body disc with a T oc R~* tem- 
perature distribution, with T = 10, 000K at the inner disc 
radius ri = 0.1-Rl and with an outer disc radius = Rl 
was used to calculate eclipse profiles in the R and B band. 
Fitting model / to these eclipses showed no significant dif- 
ference in outer disc radius between the R band and CuSC>4 
filter fits. 

The inner disc boundary and the index n in the flux 
distribution are poorly constrained. The eccentricity is ro- 
bustly non-zero; the changing eclipse shape demands a non- 
axisymmetric disc. 



The most interesting result is that all three datasets 
have 4>disc around 0.4 to 0.5. This means that in our el- 
liptical model the secondary star sweeps past the smallest 
radius part of the disc at superhump maximum - a result 
unexpected if tidal stressing of the disc by the gravitational 
influence of the secondary star is responsible for the super- 
hump light. The implications of this result are discussed 
later. 

We adjusted the model so that the eccentricity varied 
during the superhump cycle as e(4> s h) = eocos 2 Tvcf> s h. This 
will be referred to as fit e. This variation in eccentricity 
follows Simpson & Wood's (1998) simulation in which the 
disc varies between being highly eccentric at the superhump 
maximum to almost circular away from the superhump. The 
results for r m i n and a max change very little, with a max again 
larger in the red (1995) than the blue (1991 and 1993). 4>disc 
is unchanged from fit / within the errors for all three years. 
The maximum eccentricity, eo, is larger than when e was 
constant. This is expected since the eccentricity is demanded 
by the variation in O-C mid-eclipse times, and these O-C 
times are non-zero at times when e is less than eo. 

Next, we obtained fits in which the eccentricity was 
again constant, but where allowed to vary be- 

tween ai at superhump maximum to a 2 half a superhump 
period later; a max ((f> 3h ) = ai + (a 2 — ai)sin 2 ir</> sh . This will 
be referred to as fit a. This was an attempt to reproduce 
the observed variations in eclipse width (Figure []). r m i„, e 
and (f>disc are essentially the same as in the first fit, while 
the values of a\ and a 2 imply a variation in a max of ampli- 
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Table 6. Treatment of a ma ^ and e in our four models, r^in, 
and (brijvr are constant in all four models. 



Fit 




Qmax 
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constant 


constant 
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constant 


eocos 2 ir4> sh 
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h (a 2 - a 1 )sin 2 Tr4> 3h 


constant 
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ai - 


h (a 2 - a±)sin 2 ir<t> sh 


eocos 2 Tr4> sh 



February 1993 average orbital lightcurves 
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Figure 7. Average orbital lightcurves for each year. The data for 
each year are sorted into two groups, so that all eclipses in each 
group have mid-eclipse times with a similar supcrhump phase. 
The values of <f)sh for each group are shown in Table |i| The third 
curve in each plot is the difference between the two lightcurves. 



tude 20 - 25 per cent in the blue and 14 per cent in the red, 
with the disc being smallest at superhump maximum. At 
its largest, the disc extends to the edge of the Roche lobe. 
While these implied variations in disc size are large, they are 
comparable to the uncertainties in <2i and a 2 , and so must 
be treated with caution. 

The final variation on our model was to allow both e 
and a max to vary as described above. This will be referred 
to as fit b. The eccentricity and <f)disc are little different from 



30.5 31.0 

Orbital phase 



Figure 9. One night of the January 1995 data with the best 
fitting lightcurve using model a plotted as a continuous line. The 
discontinuities in the fit illustrate the different dis c flu xes and 
uneclipsed fluxes allowed for each orbit. See Section p,5[ 



the fits with e varying periodically and a max constant, while 
the values of a\ and a 2 follow those from the previous fit (e 
constant and a max varying). 

The treatment of a max and e for each fit is summarized 
in Table ^. In Figure ^ we show part of the fit to the 1995 
dataset using model a. The fit (shown as the continuous line) 
shows the different disc flux and uneclipsed flux allowed by 
our model in the form of the discontinuities at phase 30.5 
and 31.5. The level of flickering out of eclipse can also clearly 
be seen, and was taken into account in our estimate of the 
errors as described above. 

Formally, the best model is that which achieves the low- 
est value of reduced x 2 (Xr)- I n Figure [u] we show the val- 
ues of Xr achieved for each model and dataset relative to 
the lowest. The minimum Xr achieved was around 0.8 for all 
datasets and models. This figure shows that the fits a and 
6, i.e. those in which a max varies on the superhump cycle, 
produce significantly better fits to the 1993 and 1995 ob- 
servations. The variation of e during the superhump cycle 
has little effect on the quality of these fits. There is less sig- 
nificant difference between the x 2 achieved by the different 
fits to the 1991 observations, although allowing a max or e to 
vary during the superhump cycle produces a better fit than 
when they are both constant. The significant reduction in Xr 
achieved by allowing a„ lax to vary implies that this model 
best represents the behaviour of the system. 

It is also interesting to compare how well each model 
predicts the variation in the eclipse width and O-C mid- 
eclipse times. The predictions of each model are plotted in 
Figure ^. It is periodic variation in these two eclipse charac- 
teristics which requires the disc to be eccentric. The models 
poorly reproduce the O-C variations in the 1991 data. While 
the phasing of the predicted variation agrees with the ob- 
servations, the amplitude is too low. The fits in which e 
varies during the superhump cycle predict a larger modula- 
tion in O-C times, a result of the larger eccentricity in these 
fits, but the agreement for these fits is still poor. The 1991 
lightcurves suffer more from flickering than the 1993 and 
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Figure 10. The mini mum values of \ 2 achieved for the various 
models fitted (Section 3.5). 



1995 data, with many eclipses distorted as a result. This is 
the most likely explanation for the poor agreement between 
our model and the 1991 lightcurves. The agreement between 
the predicted and observed O-Cs is very good for all models 
for the 1993 observations. The variation in eclipse width is 
only reasonably modelled by those fits in which a max varies. 
The same is true of the 1995 fits. 

The result of these comparisons between the different 
models, both the formal comparison of reduced \ 2 an d the 
more subjective 'chi-by-eye' considerations of the O-C times 
and eclipse widths is that the models in which a ma x varies 
during the superhump cycle predict the observations better 
than those in which a max is constant. 

All four models agree on three important points. The 
values of a max , ai and a 2 show that the disc is large, filling 
at least about 50 per cent of the Roche lobe area. The disc 
is not axisymmetric. From the consistent values of 4>disc we 
see that when the superhump reaches maximum light, the 
light centre of the disc is on the far side of the white dwarf 
from the donor star. 



3.6 Eclipse mapping 



In Section 3.5 we used the changing eclipse profiles to con- 



strain the parameters of a model intensity distribution. An 
alternative method for investigating the distribution of emis- 
sion in the orbital plane is the commonly used eclipse map- 
ping technique developed by Home (1985). This method 
assumes that intensity distribution is fixed in the corotat- 
ing frame of the binary, lies flat in the orbital plane and 
is constant. The surface of the secondary star is described 
by its Roche potential surface. Maximum entropy methods 
are used to obtain the intensity distribution by comparing 
the calculated lightcurve and observed lightcurves. The x 2 
statistic is used to ensure consistency between the observa- 
tions and the fitted distribution, while the entropy is used 
to select the most appropriate solution from the multitude 
of possibilities. The entropy is usually defined such that the 
final solution is the smoothest or most axisymmetric map 
consistent with the observations. This technique has been 



widely used, and O'Donoghue (1990) employed it to locate 
the source of the strong normal superhumps in Z Cha. 

In order to study the shape of the precessing disc in 
V348 Pup, the PRIDA eclipse mapping code of Baptista & 
Steiner (1991) was modified so that the intensity distribu- 
tion was fixed in the precessing disc frame rather than the 
corotating frame of the bina ry. Each year's data was split 
into two groups (Section 3.4) but the lightcurves were not 



folded on orbital phase. This enabled us to obtain two maps 
for each year, corresponding to the groups given in Table ^. 
Since we expect the intensity distribution to change through- 
out the superhump cycle, grouping the eclipses as described 
means that the intensity distribution should be roughly 
the same for all eclipses in a group, an assumption of the 
eclipse mapping method. The superhump mo dul ation was 



3.5. Normal- 



subtracted from each lightcurve as in Section 
ization of the lightcurves was achieved by using the values 
for total disc flux and uneclipsed flux for each orbit obtained 
during the fitting procedure in Section [T^. The uneclipsed 
flux was subtracted from each orbit and fluxes were then 
rescaled to produce an effective disc-only lightcurve. Various 
other normalization techniques were tested, and the detail of 
the reconstructed maps was sensitive to these changes. We 
used orbital parameters q — 0.31 and i = 81°, and looked 
for the most axisymmetric solution consistent with the data. 

The most consistent result revealed by these eclipse 
maps is that the emission at superhump phase 0.5 is less 
centrally concentrated than at superhump phase 0. This 
is illustrated in Figure |ll| which shows the maps for the 
two 1995 groups of eclipses. Figures [l^ a to |l2| c show the 
azimuthally averaged brightness distribution of each map. 
They show the flux at radius r, F r , multiplied by r; this 
quantity is proportional to the total flux in an annulus at 
radius r. Figures [H^j and|l2{; show how the disc extends fur- 
ther out at superhump phase 0.5 than at phase 0, while the 
curves in Figure |l2^ , are both nearly the same, as expected 
since both of these curves show the situation roughly half 
way between superhump phase and 0.5. This result is in 
agreement with the results of our fits of model a in which 
the disc size was allowed to change. These fits showed that 
the size of the emission region is larger at superhump phase 
0.5 than at phase 0. The maps are asymmetric, but due to 
the sensitivity described above, we draw no conclusions from 
the detailed structure. 



4 DISCUSSION 

The phase of the superhump relative to the conjunction of 
the line of centres of the system and the semi-major axis 
of the disc should make it possible to determine whether 
the bright spot model or the tidal heating model better 
explains the source of the superhump. The simplest tidal 
model predicts that the superhump light should peak when 
(or slightly after) the largest radius part of the disc coincides 
with the line of centres. This is because the tidal interaction 
is strongly dependent on distance from the secondary, and 
so will be most significant in regions where the disc extends 
out close to the LI point. However, if the bright spot model 
is to be believed, then the the superhump light source will 
be brightest when the accreting material has the furthest to 
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Table 7. P arameters resulting from fitting various models to the lightcurves (see 
Section 3.5). For the behaviour of a max and e in each model see Table pi 



Fit Year 



i/Rh dmax/Rh e 

ax/R L a 2 /R L e 



/ 


1991 


0.13(13) 


0.81(10) 


0.035(03) 


0.61(81) 


0.46(17) 


/ 


1993 


0.12(09) 


0.83(07) 


0.106(02) 


0.94(65) 


0.41(11) 


/ 


1995 


0.11(13) 


0.90(09) 


0.082(03) 


0.64(71) 


0.53(13) 


e 


1991 


0.13(10) 


0.82(08) 


0.121(06) 


0.62(76) 


0.35(15) 


e 


1993 


0.09(10) 


0.81(07) 


0.144(03) 


0.67(60) 


0.41(03) 


e 


1995 


0.08(10) 


0.88(06) 


0.102(04) 


0.48(52) 


0.55(15) 


a 


1991 


0.12(14) 


0.72(12) 


0.88(17) 


0.029(02) 


0.59(84) 


0.51(22) 


a 


1993 


0.11(11) 


0.76(08) 


0.99(16) 


0.099(04) 


1.11(73) 


0.41(15) 


a 


1995 


0.10(16) 


0.78(11) 


0.91(09) 


0.068(05) 


0.34(96) 


0.53(21) 


b 


1991 


0.12(15) 


0.74(09) 


0.89(20) 


0.104(08) 


0.55(79) 


0.35(18) 


b 


1993 


0.10(09) 


0.76(08) 


0.98(13) 


0.140(05) 


0.99(61) 


0.41(12) 


b 


1995 


0.09(17) 


0.80(09) 


0.93(10) 


0.084(06) 


0.42(61) 


0.56(22) 
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Figure 12. Azimuth allv averaged flux in maximum entropy eclipse maps. Quantity plotted, rF r , is proportional to flux in annulus at 
radius r. See Section |3.6|. 
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Figure 11. Maximum entropy eclipse maps of the intensity dis- 
tribution in the precessing frame of the disc from the 1995 data. 
The average orientation o f the primary Roche lobe for each group 
is also shown. See Section 3.6. The grey scale is the same for both 
maps. 



fall. In other words, the superhump should occur when the 
stream impacts on the disc at its smallest radius. 



The mid-eclipse times shown in the top row of panels 
in Figure ^ show the eclipses to be earliest around super- 
hump phase 0.75 in all cases. Assuming that the centre of 
light of the eccentric disc is offset from the white dwarf in 
the direction of the largest radius, we can deduce the disc 
orientation during these eclipses to be as shown in Figure 
p^a. A quarter of a superhump period later, the orientation 
of the disc has barely changed, the secondary will be lined 
up with the smallest radius part of the disc and the super- 
hump phase will be 0.0 (Figure [l^b). Therefore superhump 
maximum occurs when the secondary star is lined up with 
the smallest radius part of the disc. The values of (j>disc in 
Table (?) agree with this deduction. This phasing is consistent 
with the bright spot model for the superhump emission but 
is in cons istent with the simple tidal heating model. In Sec- 
tion 3.4, we noted that the superhump light appears not to 
be emitted isotropically: the superhump is strongest when it 
occurs around orbital phase 0. This is easily explained if the 
major contribution to superhump light is the bright spot: 
the bright spot is most visible when it is on the nearside 
of the accretion disc. Schoembs (1986) observed late super- 
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Figure 13. The first row of Figure Q shows the eclipses to be 
earliest at superhump phase </> s h ~ 0.75, which implies orienta- 
tion (a), from which we deduce the relative phasing of disc and 
secondary star at superhump maximum (<ji s h = 0) shown in (b). 



humps in the eclipsing SU UMa dwarf nova OY Car, also a 
high inclination system. When a superhump was coincident 
with a pre-eclipse orbital hump, the combined amplitude 
was greater than that predicted for a linear superposition 
of the individual amplitudes i.e. OY Car's late superhumps 
were strongest around orbital phase 0. However, van der 
Woerd et al. (1988) studied the dwarf nova VW Hyi, con- 
cluding that there was no correlation between the orbital 
phase and amplitude of late superhumps. Since VW Hyi has 
an intermediate inclination, ~ 60° (Schoembs & Vogt 1981), 
the bright spot visibility need not vary with phase, so if the 
bright spot is the main superhump light source we expect 
no variation in superhump amplitude with orbital phase. 

Krzeminski & Vogt (1985) studied OY Car during a 
super-outburst and through variations in the O-C eclipse 
timings deduced the presence of an eccentric disc with phas- 
ing similar to that in V348 Pup. Krzeminski & Vogt's defini- 
tion of O-C time was criticized by Naylor et al. (1987), with 
Naylor et al. concluding that the O-C evidence was weaker 
than originally thought. 

Schoembs (1986) followed OY Car from early in a super- 
outburst almost until the return to quiescence, observing the 
~ 180° phase change from normal superhumps around the 
height of the outburst to late superhumps during the decline 
of the super-outburst. Patterson et al. (1995) observed the 
same change in superhump phase late in a super-outburst 
of V1159 Ori. Hessman et al. (1992) studied OY Car at the 
end of a super-outburst. By looking at the varying hot spot 
eclipse ingress times, and considering the trajectory of the 
accretion stream, they concluded that the disc was eccentric. 
The orientation of the disc at superhump maximum was very 
similar to that which we find in V348 Pup. 

The broad waveform of these late superhumps in OY 
Car (Hessman et al.) was also similar to the superhump 
modulation in V348 Pup. Such similarity between late su- 
perhumps in OY Car and the superhumps in V348 Pup is 
not surprising. Late superhumps in dwarf novae appear to- 
wards the end of the superoutburst, after the disc has had 
time to adjust to its high state. V348 Pup is persistently in a 



high state. Superhumps in a novalike system might resemble 
those to which the superhumps in a superoutbursting dwarf 
nova would tend if it remained in superoutburst for a long 
time. It seems likely that the mechanism responsible for late 
superhumps in SU UMa systems is the same mechanism re- 
sponsible for superhumps in V348 Pup. However, Skillman 
et al. (1998) observed strong superhumps in the nova- like 
TT Ari throughout 1997 whose waveform is triangular like 
those of normal superhumps in dwarf novae. 

There are many other studies of the disc structure in SU 
UMa stars during superoutburst. Vogt (1981) and Honey et 
al. (1988) found evidence for an eccentric precessing disc 
in Z Cha from the radial velocity variations of various ab- 
sorption and emission lines respectively. The very prominent 
normal superhump in Z Cha made it possible for Warner & 
O'Donoghue (1988) to study the location of the superhump 
light source. They found strong departures from axisym- 
metry in the superhump surface-brightness. O'Donoghue 
(1990) employed a modified eclipse mapping technique to Z 
Cha lightcurves and found the normal superhump light com- 
ing from three bright regions of the disc rim, located near the 
LI point and the leading and trailing edges of the disc, con- 
cluding that the superhumps are tidal in origin, and that 
a highly eccentric disc with a smooth brightness distribu- 
tion is not necessary to explain superhump behaviour. One 
anomalous eclipse did confine the superhump light source in 
Z Cha to the region of the quiescent bright spot, van der 
Woerd et al. (1988) concluded that the late superhumps in 
VW Hyi come from an optically thin plasma and could be 
a result of tidal interaction. 

In the SPH simulations of Murray (1996, 1998) pseudo- 
lightcurves are produced by assuming the heat produced by 
viscous dissipation to be radiated away where it is gener- 
ated. Murray (1996, 1998) reveals an extended superhump 
light source in the outer disc, while Murray (1996) reveals an 
additional superhump modulation which arises from the im- 
pact of the accretion stream with the edge of the disc occur- 
ing at a varying depth in the primary Roche potential. This 
additional weaker superhump modulation is approximately 
180° out of phase with the modulation due to tidal stressing, 
another similarity between late superhumps in dwarf novae, 
the persistent superhumps in V348 Pup and the bright spot 
model. 

If we consider the stream to impact the disc at radius r 
in a 1 potential, then the luminosity, L, of the bright spot 
should vary roughly as A(i). Considering the change in r as 
the disc with eccentricity e precesses, we get ~ 2e. The 
eccentricities we find are in the range 0.035-0.15 predicting 
superhump fractional amplitude in the range 0.07-0.3. This 
is consistent with the measured superhump amplitudes in 
V348 Pup (Table |). 

While we limit the conclusions drawn from our eclipse 
maps in Section 3.6, there are a number of points deserving 
consideration. 

Our eclipse maps do not show evidence of a bright spot, 
but this does not rule out the possibility that a bright spot is 
the source of the superhump light, for the following reasons. 
First, we subtracted the superhump modulation from the 
lightcurves before performing the eclipse mapping, which 
should reduce the contribution of the bright spot in the maps 
if it is the primary source of the superhump light. Also, 
our maps are fixed in the precessing disc frame, rather than 
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the orbital frame of the system, so the hot spot should be 
blurred azimuthally in our maps by ~ 70° corresponding 
to the eclipse width of the system of ~ 0.2 in orbital phase. 
There will be additional azimuthal blurring since the eclipses 
contributing to each map have a spread of disc orientations 
at mid-eclipse corresponding to the values of in Table 
|E| Azimuthal structure in the maps is also suppressed by 
looking for the maximally axisymmetric solution. 

The eclipse maps tell us that the azimuthally averaged 
radial extent of the emission is lowest at superhump max- 
imum, shown in Figure |l^. If this change in extent of the 
emission region is interpreted as a result of a changing disc 
size, the the smaller disc radius at superhump maximum is 
consistent with the bright spot model for the superhump 
light source. 

In the SPH models of Simpson & Wood (1998), the 
symmetry axis of the disc is aligned roughly perpendicular 
to the line of centres of the system when the superhump 
reaches maximum intensity. Inspecting their plots suggests 
that this model would lead to eclipses being earliest and 
widest at superhump phase contrary to our findings. Simp- 
son & Wood stress that their pseudo-lightcurves should be 
treated cautiously since no radiative processes were explic- 
itly considered. The difference in the mass ratio between 
V348 Pup and the values considered in Simpson & Wood's 
simulations may affect the phasing of the early eclipses and 
superhump maximum, given that the predicted superhump 
waveform is sensitive to q. Furthermore, the spiral density 
waves in their simulations complicate the structure, so that 
simulated maps of the intensity may in fact produce reason- 
able agreement with our findings. 

SPH simulations (Murray 1996 & 1998, Simpson & 
Wood 1998) show the behaviour of tidally distorted discs to 
be more complicated than a simple eccentric disc, and with 
treatment of radiative processes the predictions are likely to 
become even more complicated. Once such models are devel- 
oped further, comparisons with observation should provide 
a more complete understanding of superhump phenomena. 



ations is currently without conclusive explanation, but they 
could result from rotationally symmetrical structure in the 
disc, such as spiral waves, which have been predicted in sim- 
ulations (e.g. Simpson & Wood 1998), and directly observed 
in the dwarf nova IP Peg (Steeghs, Harlaftis and Home 
1997). They may be a result of the correlation between su- 
perhump amplitude and orbital phase. In a high inclination 
system like V348 Pup the eclipse of the superhump light 
source may also explain these signals. 

The phasing of the superhumps in V348 Pup is like that 
of late superhumps detected at the end of superoutbursts in 
SU UMa systems such as OY Car and VW Hyi (Schoembs 
1986 and van der Woerd et al. 1988). The broad form of the 
superhump and the link between superhump amplitude and 
orbital phase are also similar to late superhumps in OY Car. 
By identifying the superhump mechanism in novalikes we are 
likely to also understand the mechanism for late superhumps 
in SU UMa systems. It is possible that this mechanism could 
be different from that which produces common superhumps, 
generally accepted to be a result of tidal stresses acting on 
an eccentric disc (O'Donoghue 1990). 
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5 SUMMARY 

The eclipsing novalike cataclysmic variable V348 Pup ex- 
hibits positive superhumps. The period of these superhumps 
is in agreement with the generally observed e — P or b relation 
for superhumpers. 

Using the relation for q as a function of the superhump 
period excess, e, (Patterson 1998b) we estimate q = 0.31. 
Using the eclipse width we then estimate an orbital inclina- 
tion of i = 81? 1 ± 1?0. 

Variations in the O-C mid-eclipse times and eclipse 
widths strongly suggest that V348 Pup harbours a process- 
ing eccentric disc. We quantify this conclusion by fitting an 
eccentric disc model to the lightcurves. The relative orien- 
tation of the disc and secondary star at superhump maxi- 
mum is more easily explained by the bright spot model for 
the superhump light source than the tidal heating model. A 
correlation between the amplitude and orbital phase of the 
superhumps is also more easily explained by the bright spot 
model. 

Additional signals are detected at frequencies close to 
harmonics of the orbital frequency. The source of these vari- 
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